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CHAPTER  I 


INTRODUCTION 


1.1  Problem  Statement 

The  ability  to  analyze  error  probability  is  an  important  step  in  designing  and 
implementing  a  communication  system.  In  some  cases,  this  analysis  can  be  very  easy 
while  in  others  it  can  be  very  difficult.  There  are  numerous  techniques  available  to 
analyze  error  performance.  The  most  desirable  is  to  find  an  analytic  solution  for  the  error 
probability.  In  many  cases,  analytic  solutions  are  extremely  difficult  to  obtain  and 
hardware  implementation  or  computer-aided  simulation  may  be  required  to  facilitate 
analysis  [1],  Additionally,  numerical  and  residue  evaluation  techniques  may  be  used  in 
conjunction  with  the  analytic  and  computer-aided  methods.  The  method  of  choice  is  the 
one  that  obtains  accurate  error  performance  results  quickly,  easily,  and  economically. 

One  method  to  analyze  error  performance  of  a  communication  system  is  proposed 
by  Biglieri,  Caire,  Taricco,  and  Ventura-Traveset  in  [2]  and  [3],  The  method  makes  use 
of  the  two-sided  Laplace  transform  of  the  probability  density  function  (PDF)  of  a 
decision  statistic.  A  Gauss-Chebyshev  quadrature  rule  is  then  applied  to  the  inverse  of 
this  Laplace  transform  to  calculate  the  error  probability.  This  thesis  will  center  on  fully 
explaining  the  derivation  of  this  method.  Additionally,  the  method  is  verified  by 
analyzing  the  error  performance  of  several  communication  systems  of  varying 
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complexity.  A  major  objective  is  to  determine  and  explain  how  to  optimally  implement 
the  method  in  order  to  reduce  computational  time  while  maintaining  numerical  accuracy. 
As  shown  in  this  thesis,  the  method  has  potentially  broad  ranging  utility,  especially  when 
a  closed  form  solution  for  the  error  probability  is  not  obtainable  using  other  techniques. 

1.2  Literature  Search 

As  stated,  there  are  many  techniques  used  to  analyze  error  probability  in 
communication  systems.  These  can  be  broadly  grouped  under  finding  analytic  solutions, 
computer-aided  analysis,  and  hardware  implementation.  Furthermore,  numerical  and 
residue  techniques  may  be  used  along  with  these  broad  methods. 

Finding  an  analytic  solution  for  the  error  performance  is  the  most  desirable 
technique.  The  resulting  solution  may  have  a  closed  form  or  include  an  integral  requiring 
numerical  evaluation.  As  stated  in  [1],  analytic  analysis  is  very  useful  in  the  early  stages 
of  system  design  for  exploring  relationships  between  design  parameters  and  system 
performance.  However,  such  analysis  is  usually  based  on  simplified  models  and  analytic 
results  are  generally  harder  to  obtain  as  system  complexity  increases.  The  problem  is 
further  compounded  if  one  is  required  to  compare  the  error  performances  of  several 
systems. 

There  are  numerous  examples  of  employing  analytic  methods.  One  need  only 
look  in  a  digital  communications  book  to  find  many  such  examples.  Another  example  is 
that  of  Lindsey’s  investigation  into  error  probabilities  for  binary  andM-ary  signals  in  a 
Rician  fading  channel  in  [4],  Simon  and  Alouini  also  find  analytic  results  for  error 
performance  in  generalized  fading  channels  in  [5],  To  accomplish  this  analysis,  they 
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make  use  of  alternative  forms  of  the  Marcum  Q-function  described  in  [6]  and  [7],  This 
shows  how  the  numerical  Q-function,  which  is  an  integral  with  no  closed  form,  is 
manipulated  to  obtain  analytic  results. 

The  next  broad  category  to  evaluate  error  performance  is  computer-aided 
analysis.  This  could  range  from  simple  implementation  of  analytic  equations  to  a  full¬ 
blown  waveform  level  simulation  [1],  Computer  simulation  methods  include  Monte- 
Carlo  simulation,  importance  sampling,  extreme  value  theory  analysis,  tail  extrapolation, 
quasianalytic  techniques,  and  others.  Using  simulation  allows  analysis  of  a  system  to  any 
desired  level  of  detail.  Furthermore,  a  multitude  of  complex  systems  can  be  considered 
with  relative  ease  and  in  less  time  than  through  the  other  broad  methods. 

Several  examples  of  how  simulation  is  used  to  determine  error  performance  were 
examined.  Chen,  Lu,  Sadowsky,  and  Yao  explain  a  strategy  of  how  to  employ 
importance  sampling  in  [8],  Also,  in  [9],  Stadler  and  Roy  describe  how  to  use  adaptive 
importance  sampling  to  decrease  computational  time.  Finally,  Tranter  and  Schneider 
introduce  another  simulation  technique  called  Partial  Sum  of  Products  (ParSOP)  to 
simulate  digital  communication  systems  in  nonlinear  channels  in  [10], 

Hardware  implementation  is  another  possible  way  to  determine  the  error 
performance  of  a  system.  Basically,  this  method  involves  actually  building  a  system  and 
measuring  its  error  performance  over  the  communication  channel  being  used.  This  could 
also  simply  involve  measuring  channel  characteristics  at  the  desired  frequencies, 
bandwidths,  or  other  design  parameters  possibly  affected  by  the  channel  [1],  Two 
examples  of  channel  characterization  are  shown  in  [11]  and  [12],  Results  using  this 
method  are  often  accurate  and  reliable.  However,  the  drawback  is  that  it  is  expensive  and 
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not  very  flexible.  This  technique  is  useful  during  the  later  stages  of  the  design  process 
when  there  are  fewer  choices  to  consider. 

In  reality,  a  combination  of  the  broad  methods  of  error  performance  analysis  is 
required  when  actually  designing  a  system.  Analytic  analysis  is  the  basis  of  all  the 
methods.  This  is  because  the  theoretical  attributes  of  a  system  must  be  understood 
through  at  least  some  analytic  analysis  before  one  can  model  and  actually  implement  the 
system.  It  may  not  be  possible  to  account  for  all  of  the  actual  sources  of  errors  in  a 
system  by  analytic  analysis,  but  it  must  be  the  basis  of  any  design  process.  Once  a  basic 
understanding  of  a  system  is  in  hand,  computer  simulation  can  be  conducted  to  fully 
investigate  all  possible  sources  of  errors.  This  simulation  might  involve  using  actual 
measured  attributes  of  a  system  under  consideration  or  measured  channel  characteristics. 
Finally,  simulated  waveforms  may  be  used  in  testing  the  resulting  hardware 
implementation  and  to  help  verify  the  validity  of  the  computer  simulation  used  in  the 
design  process. 

As  previously  stated,  numerical  and  residue  integration  techniques  are  commonly 
used  in  error  performance  analysis  using  the  analytic  or  computer-aided  analysis 
methods.  Using  the  Marcum  Q-function  could  be  considered  a  “numerical  technique” 
since  it  is  an  integral  that  is  not  expressible  in  closed  form  and  must  be  numerically 
computed.  Additionally,  the  application  of  a  Gauss-Quadrature  rule  in  the  method  under 
consideration  in  this  thesis  can  be  considered  a  numerical  technique.  A  residue 
integration  technique  could  be  used  in  place  of  this  Gauss-Quadrature  rule  to  evaluate 
error  performance.  This  thesis  centers  on  using  the  Gauss-Quadrature  rule  to  aid  in 
finding  analytic  error  performance  results  and  implementing  them  using  a  computer. 
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Thus,  it  is  shown  how  a  combination  of  analytic  and  computer-aided  methods,  along  with 
a  numerical  technique,  is  used  to  evaluate  error  performance  in  communication  systems. 

1.3  Thesis  Overview 

Chapter  1  has  introduced  the  method  under  consideration  along  with  a  short 
description  of  the  broad  categories  of  error  performance  analysis  methods. 

Chapter  2  explains  the  method  derivation.  First,  the  relationship  between  the 
probability  of  error  and  the  inverse  of  the  two-sided  Laplace  transform  is  established. 
This  is  followed  by  application  of  the  Gauss-Chebyshev  quadrature  rule. 

Chapter  3  implements  the  method  to  analyze  error  performance  in  several 
communication  system  models.  Each  system  is  fully  described,  a  theoretical  probability 
of  error  is  derived,  and  the  method  is  applied.  This  is  followed  by  implementing  the 
method  using  a  computer  and  comparing  the  resulting  error  probabilities  to  the  theoretical 
results. 

Chapter  4  discusses  the  method  implementation.  Suggestions  for  implementation 
are  made  as  an  aid  to  generating  accurate  results  with  as  little  computational  time  as 
possible. 

Finally,  Chapter  5  is  the  concluding  chapter  where  recommendations  for  further 


research  are  made. 


CHAPTER  II 


METHOD  DEVELOPMENT 

This  chapter  introduces  the  method  under  consideration  that  is  used  to  obtain  error 
probability  in  communication  systems.  An  explanation  of  how  the  inverse  of  the  two- 
sided  Laplace  transform  of  the  PDF  of  a  decision  statistic  is  related  to  error  probability  is 
first  presented.  After  establishing  this  relationship,  the  application  of  the  Gauss- 
Chebyshev  quadrature  rule  to  numerically  integrate  this  inverse  transform  is  discussed. 

2.1  Probability  of  Error  using  the  Two-Sided  Laplace  Transform 

A  digital  communication  system  has  a  finite  signal  set  used  to  convey  information 
through  a  communication  channel.  At  the  receiver,  a  received  signal  is  compared  to  all 
possible  transmitted  signals  to  generate  decision  metrics.  A  decision  metric  is  basically 
the  probability  that  a  received  signal  resulted  from  one  of  the  possible  transmitted  signals. 
Decision  metrics  are  determined  for  each  signal  of  the  signal  set.  The  largest  decision 
metric  is  used  to  decide  which  signal  was  transmitted,  hopefully  with  high  probability. 

The  method  under  consideration  begins  by  defining  the  metric  m(y,x )  where  x  is 
the  signal  that  was  transmitted  and  y  is  the  signal  received.  Again,  a  decision  is  made  at 
the  receiver  by  selecting  from  all  possible  transmitted  signals  the  one  that  has  the 
maximum  metric  [2],  Assuming  that  x  was  sent,  the  pairwise  error  probability  (PEP)  is 

P{x (2.1) 
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Basically,  the  PEP  is  the  probability  that  when  x  is  transmitted  and  y  received  the  metric 
m(y,x’)  is  greater  than  m(y,x).  The  variable  x’  represents  all  other  possible  transmitted 
signals  other  than  x.  Next,  the  random  variable  7  is  defined  as 

Y  =  m(y,x)-m(y,x')  (2.2) 

where  it  is  obvious  that  the  signal  is  received  in  error  if  Y  is  negative.  Therefore,  the  PEP 
is  related  to  the  random  decision  metric  by 

P{x->x’}  =  P{7<0}.  (2.3) 

Thus,  as  shown  in  (2.3),  finding  the  error  probability  of  a  communication  system  reduces 
to  finding  the  probability  that  the  random  decision  metric  Y  is  negative.  Since  the 
decision  metric,  7,  is  a  random  variable  it  will  have  a  PDF  associated  with  it  that  defines 
how  the  variable  is  distributed.  The  method  continues  by  taking  the  two-sided  Laplace 
transform  of  this  PDF,  py(y),  resulting  in 

oo 

0)  =  E(e~sY )  =  J  e~sypY  (y)dy  (2.4) 

—oo 

where  s  =  a  +  j® . 

Based  on  Laplace-transform  theory  and  (2.4),  there  are  several  important 
observations  to  be  made  before  continuing  the  method.  First,  the  integral  of  (2.4)  has  a 
region  of  convergence  (ROC)  that  is  a  vertical  strip  in  the  complex  5-plane.  This  is 
analogous  to  the  derivation  in  [14]  concerning  analysis  of  system  stability.  The  ROC 
comes  from  the  fact  that  the  two-sided  Laplace  transform,  X(s),  of  some  function,  x(u), 
defined  by 


oo 

X(s)=  jx^e^du, 

-oo 


(2.5) 
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converges  absolutely  if 


1 


x(u)e 


<  00 


(2.6) 


where s  =  a+  ja>.  The  requirement  for  convergence  from  (2.6)  can  be  rewritten  as 


oo 

^\x{u^e~au  du  <  oo  . 

-oo 


(2.7) 


Correspondingly,  the  Laplace  transform  converges  absolutely  if  a  real  positive  number  A 
exist  such  that  for  some  real  b  and  c 


,  ,  .  f  Ae™,  u>  0 
bdwl<( 

1  V  *  \  Aebu,  u<  0 


(2.8) 


By  breaking  the  integral  (2.5)  into  the  two  parts  where  u  >  0  and  u  <  0,  it  can  be  shown 
that  the  singularities  of  X(s)  arising  from  the  portion  of  x(u)  for  u  >  0  will  lie  to  the  left  of 
the  line  s  =  c.  Similarly,  the  singularities  of  X(.s)  from  the  part  of  x(u)  for  u  <  0  will  lie  to 
the  right  of  the  line  s  =  b.  Therefore,  it  is  known  that  the  ROC  of  (2.5)  will  be  some 
vertical  strip  in  the  5-plane.  Thus,  if  the  conditions  in  (2.8)  are  met  then  the  ROC  is 
defined  by 

c  <a  <b .  (2.9) 


Again,  consider  the  Laplace  transform,  Jf(5),  of  some  function  x(u).  The 
reasoning  that  follows  is  similar  to  that  concerning  bounded-input,  bounded-output 
(BIBO)  stability  in  systems  analysis.  If 
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<  oo  (210) 

-oo 

then  any  singularities  resulting  from  the  portion  of  x(u)  for  u  >  0  must  lie  to  the  left  of  the 
jco  axis  in  the  5-plane  [13,14],  Similarly,  the  singularities  from  the  part  of  x(u )  for 
u  <  0  must  lie  to  the  right  of  the  jco  axis.  Therefore,  from  (2.9),  the  ROC  of  the  two-sided 
Laplace  transform  of  some  function  x(u)  is  a  vertical  strip  in  the  5-plane  that  includes  the 
jco  axis  if  x(u)  meets  the  requirement  of  (2. 10). 

In  the  method  under  consideration,  the  two-sided  Laplace  transform  of  the  PDF 
associated  with  the  decision  metric  Y  shown  in  (2.4)  is  desired.  From  random  variable 
theory  we  know  that  the  PDF  of  the  random  decision  metric  Y  has  the  property 

pY{y)>0  (2.11) 

for  oo  <y  <  oo.  This  leads  to  the  relationship 

\Pr  O0|  =  Py  O'),  —  00  ^  T  —  00  (212) 

and 

oo  co 

\pY(y)dy=  \\Pr(y)\dy  =  \'  (213) 

-oo  -oo 

Therefore,  (2.13)  along  with  (2.10)  imply  that  the  two-sided  Laplace  transform  of  the 
PDF  pyiy)  has  a  ROC  which  is  a  vertical  strip  in  the  5-plane  that  includes  the  imaginary 
axis.  This  vertical  strip  is  bounded  by  the  singularities  of  <J>y(5)  that  are  closest  to  the 
imaginary  axis.  The  possible  i  singularities  of  <by(5)  are  denoted  in  [2]  as  5,.  Therefore, 


the  ROC  is  ct]  <  a  <  a 2  where 
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0 

><X 

il 

e 

(2.14) 

and 

cc2  =  min  (Re[^  ]) 

Re[^  ]>0 

(2  15) 

It  is  evident  from  (2.14)  and  (2.15)  that  a\  will  have  a 

negative  value  and  a 2  a  positive 

one. 

The  next  step  in  the  method  development  involves  applying  integration  by  parts  to 

(2.4).  To  accomplish  this,  let 

n{y)  =  e~\ 

(2.16) 

u'(y)  =  du(yK-se-» 

dy 

(2.17) 

v'(y)  =  pY(y), 

(2.18) 

and 


VOO=  \p,iPW-  (2.1 9) 

—00 

The  assignment  of  (2.19)  is,  by  definition,  the  cumulative  distribution  function  (CDF)  of 
the  random  variable  Y  [15],  Therefore,  v(y)  is 

vO)  =  ^0).  P20) 

Next,  by  applying  the  integration  by  parts  formula  we  get 

(5)  =  J  u(y)v'  (y)dy  =  u(y)v(y)  -  J  v(y)u'(y)dy 

00 

=  e~syFY  OOQ  +  ^  \Fr  (y)e~sydy  (2.21) 

—00 

00 

=  e-*°FY  (cc)-esa>Fy  (-00)  +s\fy  (y)e~sydy 


—00 
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However,  random  variable  theory  says  that  Fy{ oo)  =  1  and  Fy{- o°)  =  0.  Assuming  that 
a  =  Re(.v)  >  0  results  in  e'sx  =  0.  Therefore,  the  first  term  in  the  last  line  of  (2.21)  goes  to 

zero.  Finally,  assume  that  as  y  — »  -oo,  Fyiy)  goes  to  zero  faster  than  esy  goes  to  oo. 
Hence,  the  middle  term  on  the  last  line  of  (2.21)  goes  to  zero.  The  final  result  is 

oo 

<D  Y(s)  =  s\FY(y)e-sycfy.  (2.22) 

—oo 

The  inverse  Laplace  transform  is  given  by 

^  C+  jco 

x(t)  =  —  \X(s)es‘dt.  (2.23) 

2jUc-j CO 

Applying  the  inverse  Laplace  transform  to  (2.22)  results  in 

C+ joo 

Fy(y)  =  —  1  <D  y{s)esy—  o  <c<«2.  (2.24) 

2^  Jjx  s 

Applying  the  fact  that  P{Y<0}  =  Fy( 0)  yields 

1  c+ 700  J  i  C+J° 0  J 

P{Y  <0}  =  ——;  f  (sy°  -  =  ^  (2.25) 

The  error  probability  of  (2.25)  can  be  evaluated  directly  through  residues.  However,  in 
many  cases  this  may  be  very  tedious.  Therefore,  in  this  thesis  we  apply  the  Gauss- 
Chebyshev  quadrature  rule  to  numerically  integrate  (2.25)  and  obtain  the  error 
probability. 

2.2  Applying  the  Gauss-Quadrature  Rule 

As  previously  mentioned,  this  paper  centers  on  evaluation  of  (2.25)  through  use  of 
a  Gauss-quadrature  rule.  To  apply  this  rule,  we  first  note  that  the  value  of  P{Y<0}  should 
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be  real;  therefore,  the  real  part  of  (2.25)  is  sought  [2],  The  real  part  of  a  complex  number 
after  complex  division  is  dictated  by 


Re 


a  +  bi 
c  +  di 


a  +  bi  c-  di 
c  +  di  c-  di 


ac  +  bd 
'cr+I2 


+  i 


be -ad 
c2  +d2 , 


ac  +  bd 
~cr+d2' 


(2.26) 


Substituting  s  =  c+jco  and  ds  =jdco,  and  applying  (2.26)  to  (2.25)  results  in  the  real  part 
ofP{7<0}  given  by 


j  r<Mc  +  » 
2%  [  c  +  ja 


/>{F<0}  =Re^fWv~  '  J-^do 


=  Re<  J-  ) Re ^  (C  +  (C  +  Ja,)]c-J%" i  (2.27) 


2  n 


C  +  JCO 


c-jco 


1  urc^<s>r(c  +  j(Q)]  +  (Q\m[<S>T(c  +  ja))]  d 
2n\  c2+co 2 


The  new  upper  and  lower  limits  of  integration  are 


U  = 


s-c 


s=c+jco 


c  +  JO O  —  C  _  J  GO 

j  J 


=  00 


and 


Z,= 


s=c-jo o 


c  -  joo  -  c  _  -  joo 
J  j 


=  —oo . 


Now  the  change  of  variables 


with  dco  given  by 


do  = 


(2.28) 


(2.29) 


(2.30) 


(2.31) 


is  applied  to  (2.27)  to  yield 
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nr<o}  =  i-J 

=— ] 

2*1 

=  _±J|Re{,|. 


cRefe}  +  a )  Im{<4 


c2  +(0 2 


-cdx 


J  x2  Vl  ~x2 


c  Re{g  }  +  C-~^- — —  Im  {q} 


-cdx 


Imfe} 


_  yjl-x2  (x 
dx 


V+c2-c2 


*2) 


(2.32) 


a/i^x2 


where  q  =  Oy(c  +  y'cy).  To  help  explain  the  upper  and  lower  limits  for  the  integral  of 
(2.32),  the  variable  transformation  of  (2.30)  is  plotted  in  Figure  2.1.  As  shown  in  (2.30), 
there  is  a  discontinuity  at  x  —  0.  Therefore,  disregarding  the  integrand  for  now,  the 
integral  is  broken  into  the  two  parts 

oo  0-  oo 

fda)=jdco  +  fda>.  (2.33) 

—oo  —oo  0+ 


As  seen  in  Figure  2.1,  the  range  of  co  over  [-oo.  O')  corresponds  to  the  range  on  x  being 
(O',  -1).  Similarly,  the  a>  range  (0+,co]  is  commensurate  with  the  range  on  x  of  (0+,  1). 
Therefore,  after  the  change  of  variables,  the  resulting  integration  is 


w  —i 

)dco  =  \ 


0 -x2yl\- 


0+ 

r  dx+  J- 


—  c 


i  x2yj\ -x2 


■dx. 


(2.34) 


Applying  the  fact  that 


o  a 

|/(r)A  = 


(2.35) 


to  (2.34)  results  in 
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00 


-co 


-1 


Vi-*2 

c 

yjl-x2 


dx 

dx . 


(2.36) 


The  two  separate  integrals  are  combined  to  yield  the  final  result  because  evaluation  of  the 
integrals  as  the  limit  approaches  zero  from  the  left  and  from  the  right  cancel  point-by- 
point.  Finally,  substituting  U  =  1,  L  =  -1,  q  =  Oftc  +  ja>),  and  applying  the  results  of 
(2.36)  to  (2.32)  yields  the  error  probability  in  the  form  of 

P{F<0}  =  -J-j 

2k  J. 


Res<t\ 


c  +  jc- 


.  VT^ 


Im<D, 


c  +  jc- 


dx 


.  (2.37) 


15 


1000 

800 

600 

400 

200 

< 

ffi  o 

o 

-200 

-400 

-600 

-800 

-1000 

-1  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 

x 

Figure  2.1.  Variable  Transformation  co  =  cyfl-^x* / x 
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The  final  step  in  the  derivation  involves  applying  the  Gauss-Chebyshev 
quadrature  rule  defined  on  page  889  of  [16]  as 

P^Ldx  =  J(xt)+Rm  (2.38) 

h-x 2  k= i 

where  xk  =  cos(2k-l)n/2m,  wk  =  jj/m,  and  Rm  is  the  remainder  or  error  term.  In  (2.38),  the 
integral  is  approximated  by  the  summation  of  the  weighted  values  of^x)  evaluated  at  the 
abscissas  xk,  where  k  =  1  to  m,  added  to  the  error  term  [20].  Applying  the  Gauss- 
Chebyshev  rule  of  (2.38)  to  the  two-sided  Laplace  transform  in  the  integral  of  (2.37) 


results  in 
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Vl-x  Jl- cos  [(2A: - \)n  1 2m\ 

Oy  c  +  jc -  =  <E>r  c  +  jc- - — v - =j — 

x  J  co^fclk -\)x  1 2m\  . 

=  f  *Jsm2[(2k-\)7t/2m])  f  sin[(2A:  -  l>r/ 2/^]^ 

cos[(2A:  -  l);r  /  2/w]  cos[(2&  -  l);r  /  2m]J 

=  Oy  (c  +  y'c  tan[(2&  - l);r  /  2/w])  =  <f>y  (c  +  jczk ) 
where  %  =  tan[(2£-l)7z/2m]  The  square  root  of  sin2(x)  equals  sin(x)  as  long  as  x  is 
restricted  to  the  range  (0,  it).  This  requirement  is  met  because  of  the  way  k  and  m  are 
defined.  Finally,  substituting  (2.39)  along  with  m  =  n!2  and  w,  =  rim  into  (2.38)  results 


P{r  <  0} = 6 + jc*>  )] 

=  - ReI°r  (c  +  JCTk)]+Tk  lm[07 (c  +  jczk)] }  +  En 


(2.40) 


where  E„  is  the  error  term.  Also,  n  is  the  number  of  nodes  and  is  assumed  to  be  even. 

A  disparity  between  the  development  above  and  the  development  in  [2]  concerns 
the  Tic  term.  In  the  final  result  of  [2],  r*  has  2 n  in  the  denominator,  while  in  the 
development  here  there  is  an  n  term  in  the  denominator.  The  derivation  given  here 
appears  to  be  correct  since,  as  n  approaches  infinity,  the  term  zu  approaches  the  function 


This  is  evidenced  in  Figures  2.2  and  2.3.  Note  that,  as  the  number  of  nodes  increases,  the 
correct  Tk  approaches  the  form  of  the  function  shown  in  (2.41). 
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The  error  term  of  (2.40),  En,  decreases  as  the  number  of  nodes  gets  larger  and 
should  disappear  when  n  =  oo.  The  variable  c  determines  the  number  of  nodes  required  to 
obtain  a  desired  degree  of  accuracy.  It  is  stated  in  [2]  that  through  applying  the  method,  a 
good  value  of  c  is  that  value  which  minimizes  <J>y(c).  Another  way  of  choosing  c,  as 
stated  in  [3],  is  half  the  value  of  the  positive  singularity  of  closest  to  the  imaginary 
axis.  One  of  the  major  objectives  of  this  thesis  is  to  determine  exactly  what  makes  these 
choices  of  c  good.  Basically,  a  good  choice  of  c  would  enable  accurate  numerical  results 
with  the  fewest  nodes  possible.  This,  in  turn,  would  result  in  less  computational  time. 
With  the  method  suitably  defined,  we  proceed  to  test  it  on  several  communication 
models. 


CHAPTER  IH 


ERROR  PROBABILITY  ANALYSIS  IN  COMMUNICATION 
SYSTEMS  USING  THE  GAUSS-QUADRATURE  METHOD 

The  method  presented  in  Chapter  II  is  now  used  to  analyze  the  error  probability  in 
several  communication  models  of  varying  complexity.  This  analysis  is  accomplished  to 
verify  the  utility  of  the  method  for  different  models.  Antipodal  baseband  signaling  in 
additive  white  Gaussian  noise  (AWGN),  BPSK  in  AWGN  with  phase  error,  BPSK  in 
AWGN  and  Rayleigh  fading,  and  non-coherent  frequency-shift  keying  (NFSK)  in 
AWGN  and  Rayleigh  fading  are  all  considered  in  this  chapter. 

3.1  Antipodal  Baseband  Signaling  in  AWGN 

3.1.1.  Introduction.  The  first  model  considered  is  the  simple  case  of  antipodal 
baseband  signaling  in  AWGN.  A  description  of  the  model  is  first  presented  followed  by 
a  derivation  of  the  theoretical  error  probability.  Next,  the  two-sided  Laplace  transform  is 
found  along  with  the  error  probability  using  (2.40)  and  the  results  are  compared  to  the 
theoretical  error  probability. 

3.1.2.  System  Description.  The  antipodal  baseband  signaling  set  consists  of  the  two 
signals 


|  0<t<Tb 


sx  ( t )  =  +A 
s,(l)  =  -A 


(3.1) 
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Essentially,  si(t)  is  sent  for  the  bit  of  1  and  S2(t)  is  sent  for  the  bit  of  0.  As  shown,  each 
signal  lasts  for  a  duration  of  Tb  seconds.  Here,  it  is  assumed  that  the  signals  occur  with 
equal  probability  although  they  need  not.  The  transmission  channel  is  influenced  by 
AWGN  that  has  zero  mean  and  a  two-sided  power  spectral  density  (PSD)  ofJVo/2.  The 
receiver  used  to  detect  the  transmitted  signal  is  depicted  in  Figure  3.1. 


A+n(t) 

or  - > 

-A+tt(t) 


t  =  n 

^  F  > 

Decision 

THRESHOLD  =  0 

Figure  3.1  Antipodal  Baseband  Signaling  Receiver 


The  output  of  the  integrator,  Y,  is 


Tb  Tb  Tb 

Y  =  f(±A  +  n(t))dt  =  j(±A)dt  +  jn(t)dt  =  ±ATb  +N.  (3.2) 


For  this  system,  the  random  variable  input  to  the  threshold  detector  has  a  Gaussian 
distribution  because  the  additive  noise  is  assumed  to  be  Gaussian.  In  general,  the 
probability  density  function  (PDF)  of  a  Gaussian  random  variable,  Z,  with  mean  mz  and 
variance  oi  is 


Pz  0)  = 


i 

PYfi 

-{z~mz)2 

r . 7  exP 

y2  naz 

2o2z 

(3.3) 


Since  n(t )  is  Gaussian  with  zero  mean  and  PSD  No/2,  the  mean,  m^,  and  variance,  <jn2,  of 
the  noise  component,  N,  from  (3.2),  regardless  of  the  signal  sent,  are 
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mN  =  E{N}  =  7-!  j  n(t)dt  l  =  J  E{n(t)}dt  =  0 . 


(3.4) 


a2N  =£{v2}-£^=£|||«(/K^)^,'*J 


Tb  Tb  \T  Tb 

=  J*  j*  jE{/i(f)w(/,)}d5fV3&  =  — -  jdt  = 


AA '? .  w0r6  • 


(3.5) 


0  0  A/ 


The  random  variable,  7,  is  conditionally  Gaussian  with  mean  and  variance,  conditioned 
on  the  signal  sent,  of 

mr\±A  =Ety\±A}=E{±ATi+N)  =  E{±ATb}+E{N}  =  ±ATb  (3.6) 


=±ATh  =0 


a 


Y\±A 


=  e^2\±a}-e2^\±a] 

=  e\a2T2  ± 2 ATbN  +  N 2 }-  A2T2 
=  e\a2T2  }±  {2ATbN}+ £{a^]-  A2T2 


=A1T} 


KoT„/2 


=  2ATiElN}+^E=^YL 


(3.7) 


=0 


The  probability  of  error  can  now  be  found  since  the  PDFs  of  7  and  TV  are  fully  described. 


3.1.3  Probability  of  Error.  At  the  receiver,  a  decision  that  +A  was  sent  is  made  if  the 
sampled  integrator  output  at  t  =  Tb  is  greater  than  0,  and  the  decision  that  -A  was  sent 
occurs  if  the  integrator  output  is  less  than  0.  With  the  probability  of  each  bit  occurring 
half  of  the  time,  the  total  probability  of  error,  Pe,  is 

PE  =  ip(7  <  0 1 +A  sent)+-^P(7  >  0 1  -,4  sent).  (3.8) 

First,  the  probability  of  error  given  that  +A  was  sent  is  found  to  be 
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P(Y< 0 1 +A)  =  P((A Tb+N)<0\+A  sent)  =  P(N  <-ATb) 

By  Symmetry 


-ATh 


ATh 


=  \pN(n)dn  =  \pN{n)dn 

—00 

1  00 

=  .  t  [  exp 


(3.9) 


n 


Aifjji). 


dn 


Next,  a  change  of  variables  of  u  =  n/(NoTb/2)u 2  is  performed  on  (3.9)  to  result  in  an  error 
probability  of 


1  °° 

P(Y  <0\+A)=-f=  f  exp 
V2tt  — — 


2  A% 
Wn 


—  u 


\du  =  Q 


2A2L 


Nn 


(3.10) 


2Eb  ) 

f  PO 

=  0.5erfc 

IV".  J 

Ilk  J 

=  0 

where  the  Q-function,  Q(u),  and  complementary  error  function,  erfc(w),  are 

Q(u)=- j=]e~‘2,2dt 
<2 n  Ju 

2  00 

erfc(«)=  -7=  jV'2 dt . 
dx  I 

The  Q-function  and  complementary  error  function  are  related  by 


(3.11) 


(3.12) 


0(M)  =  ^erfc 


(  /  A 

u. 


2) 


(3.13) 


Using  a  similar  derivation,  the  error  probability  when  -A  is  sent  can  be  found  to  be  the 
same  as  that  when  +A  was  sent.  Therefore,  the  total  probability  of  error  is 


pe=Q 


2Ef 


"A 


V  V  ^0  j 


=  0.5erfc 


\\No  ) 


(3.14) 


With  a  theoretical  probability  of  error  in  hand,  the  Gauss-Quadrature  rule  is  now  applied. 
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3.1.4  Application  of  the  Gauss-Chebyshev  Quadrature  Rule.  To  apply  the  Gauss- 
Quadrature  rule;  we  need  the  two-sided  Laplace  transform  of  the  PDF  associated  with  the 
decision  variable,  Y,  given  +A  was  sent.  The  two-sided  Laplace  transform  of  a  Gaussian 
random  variable  with  mean  zwyand  variance  ay  is  given  by 


<M^)  =  \pjWv<ty 


yj^ra 

__1 _ 

■^Ina2 

1 


2 

Y  —oo 


yJlTta 


oo 

jexp 

-00 

00 

jexp 

-oo 

00 

{exp 


2 

Y  —oo 


exp[-5y]^ 


-iy~mr)2 
2(7  y 

-y2  -m2  +  2mYy -2ajsy 
2a 2 

y 2  +  {2m r  -  2aYs)y  +  m 
2al 


dy 


2  Y 

7  ]\dy 

J 


In  [18],  the  following  definite  integral  is  given: 

{ exp[-  {ax2  +bx  +  c)]a6c  =  J—  exp[(&2  -  4 acjf  4 a  . 

-oo  ^  ^ 

Now,  applying  this  integral  to  (3.15)  results  in 


1  1 

n 

j|  2 no]  \ 

(1/2 <r}) 

4 mY  +  4cr yS  -8 mY<jYs  4 mY 


4at 


4  °Ayj 


r  1  ' 


. 

[2a2  J 

=  exp 

4a yS2  -  Smya2s  2 a2 

=  exp 

0.5 aYs2  -mYaYs 

4aY 

4  J 

- 1 

b 

=  exp 

0.5 aYs2  -/«rs] 

(3.15) 


(3.16) 


(3.17) 


For  antipodal  baseband  signaling,  the  mean  and  variance  given  that  +A  was  sent  were 
shown  in  (3.6)  and  (3.7)  to  be  mY  =  ATb  and  ay2  =  NoTb/2,  respectively.  Therefore,  the 
two-sided  Laplace  transform  for  this  case  is 

O,  (.V)  =  exp(o  25NJ,/  -ATts) 


(3.18) 
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The  Gauss-Chebyshev  quadrature  evaluation  of  the  error  probability  for  antipodal 
baseband  signaling  was  implemented  using  MATLAB  and  the  two-sided  Laplace 
transform  given  by  (3.18).  The  results  obtained  are  compared  in  Figure  3.2  with  the 
theoretical  probability  of  error  from  (3. 14).  As  shown  in  Figure  3.2,  the  theoretical  results 
obtained  using  the  Gauss-Quadrature  method  closely  match  those  computed  using  the 
closed  form.  Furthermore,  only  50  nodes  were  used  to  obtain  results  with  a  maximum 
percent  difference  of  0.13%.  Finally,  the  code  used  to  generate  these  results  is  given  in 
Appendix  A. 


Figure  3.2  Error  Probability  for  Antipodal  Baseband  Signaling 
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3.2  BPSK  in  AWGN  and  Imperfect  Phase  Reference 

3.2.1  Introduction.  This  section  considers  BPSK  in  AWGN  and  an  imperfect  carrier 
reference  in  the  receiver.  The  phase  error  resulting  from  the  imperfect  carrier  reference  is 
assumed  to  have  a  Tikinov  distribution  [17],  Once  again,  the  system  model  is  presented 
followed  by  a  derivation  of  the  error  probability.  The  Gauss-Chebyshev  quadrature  rule 
is  then  applied  to  analyze  error  performance. 


3.2.2  System  Description.  The  signal  set  for  the  BPSK  model  is 

s ,  (t)  =  -A.  cos(ry  j)l 
Y\  ;  o <t<Tb. 

s2{t)=  Ac  cosfrvjj 


(3.19) 


Again,  it  is  assumed  that  the  two  symbols  have  an  equal  probability  of  occurrence  and  the 
AWGN  has  a  mean  of  zero  and  two-sided  PSD  Nq/2.  Additionally,  a  coherent  reference 
is  required  at  the  receiver  and  it  is  assumed  that  the  phase  tracking  loop  is  imperfect; 
therefore,  there  is  a  phase  error,  (f>,  in  the  local  reference.  The  receiver  is  depicted  in 
Figure  3.3. 


2c,os{coct  +  ^) 


Figure  3.3  BPSK  Receiver  with  Phase  Error  <j) 

The  output  of  the  integrator  of  Figure  3.3  is 
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1b 

=  1 2  cos (coct  +  (ff+  Ac  cos  coct  +  n(t))dt 

0 

Tb  h 

=  |2  cos(<yc/  +  Ac  cos (oci)dt  + J 2 n(t) cos (coct  +  <f)dt 


(3.20) 


±  Ac  cos{lcoct  +  <f)±Ac  cos^ 


N 


dt  +  N  =  ±  AcTb  cos<f>  +  N 


Assuming  (j>  is  known,  Y  is  again  a  conditionally  Gaussian  random  variable.  The  mean 
and  variance  of  the  Gaussian  random  variable  N,  given  <j),  are  found  to  be 


m 


m 


f h  ]  h 

E{N }  =  jEx  J 2  cos(<yct  +  f)n(t)dt  l  =  J  2  cos (coct  +  <f)E{n(t)}dt  =0  (3.21) 

U  J  0  ^ 


and 


a 


m 


e\n2}-e2{n) 


ThTh 


=  4  cos  (coj  +  <f>)cos{(Dj'+<fi)n(t)n(f)dt'dt  > 

.0  0  J 


ThT. 


=  J 1 4  cos {a>ct  +  f)cos{coc  t'+(f)E{n{t)n(t'  )}dt'dt 


(3.22) 


0  0 


4  Nn 


JE-sv-c) 


Jb 

|  cos2  (coct  +  <f)dt 


AN, 


0.5  +  0.5  cos(2<yct  +  2<f) 


AN  T 
dt  =  =  N  T 


0 


For  simplicity,  let  -1  represent  Si(t)  and  +1  represent  S2(t).  With  the  parameters  of  N 
known,  the  mean  and  variance  of  Y,  given  <f>,  are 


m 


r \<fi. 


±,  =E{Y)  =  E{±AcTt  cos* +  N) 

=  E{±AcTt  cos^}+£{iV}  =  ±AcTb  cos^ 


(3.23) 


and 
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<x^±1  =  ^7  -  E{Y  f  }=  E\±  AcTb  cos </>  +  N  +  AcTb  cos^)2 

=  4\r2}=jv; 


(3.24) 


3.2.3  Probability  of  Error.  It  is  again  desired  to  find  the  total  probability  of  error  as  in 
(3.8).  For  this  case,  the  phase  error,  </>,  is  assumed  to  be  a  given.  The  conditional  error 
probability  is 

V*  =^(J'>oU,(O,4>  =  f»)+|f(J'<0U!(O,<6  =  «t).  (3.25) 


Assuming  that  the  signal  si(t)  was  sent,  the  probability  that  s2(t)  is  received  is  found  to  be 

dn 


_  1  2 
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=  —  erfc 
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=  —erfc 
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V27r<7w  ^ 

Jexp[- =  —erfc 


00 

r 

1 

to 

_ 1 

J  exp 

\  COS  <f> 

.  2  dl_ 

AcTb  cos  <f> 
J2o\j 


AcTb  cos</> 


2N0Tb 


COS  (j) 


=  —  erfcl 
2 


\A2X 

-  cos  <b 
'  2  Nn 


Nn 


COS0 


(3.26) 


The  substitution  Ac2Tb/2  =  Eb  was  made  in  (3.26).  The  probability  of  error  when  52(0  is 
sent  can  be  shown  to  be  the  same  as  that  of  (3.26).  Therefore,  the  total  probability  of 
error  conditioned  on  ^  is 


PE\*=t  =  °-5erfcl 


Eh 

\2Eh 

-Ecos0 

=  Q 

'  cos  (j) 

Lv^o 

_  V  ^0 

(3.27) 


The  unconditional  error  probability  is  found  according  to  [15]  by  multiplying  the 
conditional  error  probability  by  the  marginal  PDF  of  O  and  integrating  over  <j)\ 
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\PE\^ 

-00 

The  phase  error  is  assumed  to  have  a  Tikinov  PDF  defined  in  [17]  as 


(3.28) 


<  n 


(3.29) 


where  h(p)  is  the  modified  Bessel  function  of  order  zero.  The  variable  p  is  the  tracking 
loop  SNR  defined  to  be  equal  to  the  inverse  of  the  phase  error  variance,  cr/.  Finally,  the 
total  probability  of  error  is 

exp(/3  cos^) 


Pe  =  fpE\<t,=tP*(</>W=  Jrerfc 


r  1 

Eu 

f-erfc 

— b-  COS  (j) 

J  2 

-n 

_ V  -^o 

27d0(p) 


-d<j> 


— Vllerfc  aP 

4^o  (p)1  Wo 


cos^ 


(3.30) 


exp(pcos^)t/^ 


A  closed  form  solution  for  (3.30)  is  not  possible.  Therefore,  it  must  be  numerically 
integrated. 


3.2.4  Application  of  the  Gauss-Chebyshev  Quadrature  Rule.  It  is  again  assumed  that 
the  phase  error  is  a  given  and  now  the  conditional  two-sided  Laplace  transform  is  found. 
But,  since  the  conditional  PDF  of  Y  when  given  <p  is  Gaussian  with  mean  and  variance 
given  in  (3.23)  and  (3.24)  respectively,  the  two-sided  Laplace  transform  is  known  from 
the  antipodal  case  to  be 


®7i$=^(,s) —  —  ^]—  exp(o.5(77|^s  m7„s). 


(3.31) 


To  find  the  unconditional  Laplace  transform,  the  relationship  given  in  [15]  is  utilized: 

oo 

£fe(X,f)}=  E{E\g(x,r]x  =  x]}=  jE^ixjjx  =  xjpj.  (*)*  (3.32) 
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Assuming  s2(t)  is  sent  and  a  Tikinov  phase  error  distribution,  the  relationship  in  (3.32)  is 
used  to  find  that  the  unconditional  Laplace  transform  is 


=4""')=  J4""1 

—00 

=  }exp[o.5 a2r^s2  -mr{ 

=  eXP^~  w°y  ^  J exp[cos^(-  AXv  +  p)]cty 


2^4  W 

The  following  relationship  is  defined  in  [19]: 


(3.33) 


7L 

| exp[y0  COS (f\l<f>  -  2 7vl0  (/?) . 


(3.34) 


Applying  this  relationship  to  (3.33)  allows  the  two-sided  Laplace  transform  to  be 
expressed  in  closed  form  as 

a,  w  =  r.(p-^.r,)«pb-5^l.  p.ss) 

hiP) 


MATLAB  is  again  used  to  calculate  error  probabilities.  Several  values  of  loop 
SNR  are  considered.  The  loop  SNR  is  set  at  a  desired  constant  level  above  Eb/No.  For 
example,  the  PLL  SNR  is  maintained  at  0  dB  above  Eb/No  in  the  graph  of  Figure  3.4.  The 
results  are  depicted  in  Figures  3. 4-3. 6.  Time  savings  between  numerical  integration  with 
MATLAB ’s  quad8  command  and  the  method  under  consideration  are  noted  at  the  top  of 
each  graph.  Additionally,  each  graph  also  has  a  plot  of  the  error  probability  if  the  phase 
error  were  zero.  As  shown,  the  results  for  both  methods  match  to  within  1%.  More 
importantly,  there  is  up  to  99%  computational  time  savings  realized  when  using  the 
Gauss-Chebyshev  quadrature  method  as  compared  with  using  the  quad8  command.  The 
code  used  to  create  these  results  is  shown  in  Section  2  of  Appendix  A. 
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Nodes:  280;  PLLSNR:  1 0  dB  above  E /N  ;  Time  Savings  99%  ;  Max.  Percent  Difference:  0.75% 

b  o 


Figure  3.6  Error  Probability  for  BPSK  with  Loop  SNR  of  10  dB 

Above  Eb/No 

3.3  BPSK  in  AWGN  and  Rayleigh  Fading  Channel 

3.3.1  Introduction.  The  next  model  under  consideration  is  BPSK  in  a  channel 
influenced  by  AWGN  and  Rayleigh  fading.  The  Rayleigh  fading  is  assumed  to  be 
frequency  non-selective  and  to  vary  slowly  over  a  bit  period.  Again  the  system  model  is 
described,  the  error  probability  is  derived,  and  the  Gauss-Quadrature  rule  is  applied. 

3.3.2  System  Description.  The  signal  set  for  this  model  is 

®i(0  =  -^cos(<ofl0l  n<f<T 

s2(0  =  Ac  cos(«v)J 


(3.36) 
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The  signal  set  is  obviously  the  same  as  that  in  the  other  BPSK  case  and  is  reprinted  here 
for  convenience.  Again,  it  is  assumed  that  the  two  symbols  have  equal  a  priori 
probability  and  the  AWGN  has  a  mean  of  zero  and  two-sided  PSD  NJ2.  The  Rayleigh 
fading  in  the  communication  channel  is  a  random,  multiplicative  process  that  causes  the 
received  signal  to  have  a  random  amplitude.  This  process  is  assumed  to  be  frequency 
non-selective  so  it  affects  all  frequencies  in  the  same  manner.  Additionally,  the 
assumption  that  the  Rayleigh  fading  is  slow  results  in  the  received  signal  being  constant 
over  a  bit  interval.  This  allows  the  phase  to  be  estimated  without  error  [21],  Figure  3.7 
shows  the  receiver  used  in  this  model. 


The  variable  r  is  the  Rayleigh  random  variable.  In  a  similar  manner  as  that  of  the 
previous  BPSK  case,  the  output  of  the  integrator  is  found  to  be 

Y  =  ±rAcTb+N.  (3.37) 

The  Gaussian  random  variable,  N,  has  a  mean  and  variance  of 


[h 


mN  =  E{N)  =  e\  ^2cos{co ct)n(t)dt  1  =  |2 cos{o)ct)E{n(t)\dt  =  0  (3.38) 


lo 


=0 


and 
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TtTb 

<Tj=£{w2}-^M=jf4  cos(coj)co$  (ojJ')E  {n(i)n(t'  )}dt'dt 


=o 


0  0 


(3.39) 


T 

=  I  COS2  (aj)dt  =  N„Tb 

Again,  si(t)  is  represented  by  -1  and  s2(0  by  +1.  The  mean  and  variance  of  Y, 
conditioned  on  r,  are 


=E{Y}  =  E{±rAJb+N}=±rAX  (3.40) 


and 


=  e\y  -  E(Y)f\=  E\±rAcTt  +NTrAcTlf]=  N0Tb .  (341) 


3.3.3  Probability  of  Error.  In  a  similar  manner  as  in  the  previous  case,  the  conditional 
probability  of  error  is 


0.5erfc 


0.5erfc 


SAX 

2  N0 


(3.42) 


where  yb  =  Eb/No=  0.5r2Ac2Tb/No  is  the  SNR.  To  find  the  unconditional  error  probability, 
the  product  of  the  conditional  error  probability  and  the  marginal  PDF  of  R2  is  integrated 
over  r2.  Since  r  is  Rayleigh  distributed,  r2  has  a  chi-square  distribution  with  two  degrees 
of  freedom  [21],  Similarly,  yb  also  has  a  PDF  which  is  chi-square  with  two  degrees  of 
freedom  and  is  represented  by 

Pi/b  )=~\~e  b//y'b  >  Yb  >  0  (3-43) 

Yb 

where  y’b  =  ybEir 2)  is  the  average  SNR  and  Eir2)  is  the  average  value  of  r2.  Therefore, 
the  unconditional  error  probability  is 
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r  "V. 
[e  /r‘ 

~P=  ?<?  “2 12 du 

J 

0 

dyb  .  (3.44) 

ly2*  fa  J 

With  a  variable  substitution  of  z  =  ul(2yb)m  and  interchanging  the  order  of  integration, 
the  error  probability  is 
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Applying  the  integral 
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from  [18]  to  (3.44)  yields 
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r(o.5+i) 
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Finally,  applying  the  indefinite  integral 

dx 


1; 


2{ax  +  b ) 


(ax2  +bx  +  c$'2  (4  ac-Z>2)>/ax2  +  &x  +  c 

from  [18]  to  (3.47)  results  in  the  error  probability 
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PE  = 


2  Y\ 
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'l  +  A 


(3.45) 


(3.46) 


(3.47) 


(3.48) 


(3.49) 


3.3.4  Application  of  the  Gauss-Chebyshev  Quadrature  Rule.  The  decision  variable, 
Y,  of  (3.37),  given  the  variable  r  and  that  s2(0  was  sent,  has  a  Gaussian  distribution  with 
mY  =  rAcTb  and  ay2  =  NoTb.  Therefore,  the  conditional  two-sided  Laplace  transform  has 


the  familiar  form 


o„r  (.v) = expfsAy/;.!-'  -  ?;•.-) 


(3.50) 
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The  PDF  of  r  is  given  in  [21]  to  be 


pR(r) 


2  r 


^)exTW)l 


r  >0. 


(3.51) 


The  unconditional  Laplace  transform  is  then 

,  /  x  r  /  \  /  s  ,  2exp[o.5A^or,52]f 

<M4  =  f  *„,(»)»»*■=  PL^*r 


rexP|  ~TTT\r2  -AoTbsr 


l~WY 


dr.  (3.52) 


Applying  the  integral 
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(3.53) 


from  Appendix  B  to  (3.52)  yields 


3>y  (5)  =  e < 


l-e 


EJhS1 


w  2 


-EbTbs2  Jexp^  2  '  dx 


(3.54) 


where  E'b  =  E{r2^Eb  is  the  average  bit  energy.  The  last  term  of  (3.54)  could  not  be 

found  in  closed  form  and  will  be  computed  numerically. 

The  resulting  Laplace  transform  of  (3.54)  is  now  implemented  in  MATLAB  using 
the  Gauss-Chebyshev  quadrature  method  to  obtain  the  error  probability.  This  is 
compared  to  the  error  probability  from  (3.49)  in  Figure  3.8.  Figure  3.8  shows  that  less 
than  a  one  percent  difference  is  obtained  with  30  nodes.  However,  the  method  took 
considerably  longer  to  implement  than  the  closed  form  solution  due  to  the  numerical 
integration  required.  The  transform  of  (3.54)  can  be  expressed  in  terms  of  a 
complementary  error  function  and  this  would  likely  decrease  computation  time. 
However,  the  argument  of  the  complementary  error  function  is  complex  and  MATLAB 
can  not  implement  this.  Additionally,  a  method  to  find  a  good  value  for  the  constant  c 
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was  successful.  This  will  be  discussed  further  in  Chapter  IV.  The  code  again  is  shown  in 
Appendix  A,  Section  3. 


Nodes:  30;  Max .  Percent  Difference:  0.95% 


Figure  3.8.  Error  Probability  for  BPSK  in 
AWGN  and  Rayleigh  Fading 


3.4  NFSK  IN  AWGN  AND  RAYLEIGH  FADING  CHANNEL 

3.4.1  Introduction.  The  final  model  tested  is  NFSK  in  AWGN  and  Rayleigh  fading. 
Again,  the  Rayleigh  fading  is  assumed  to  be  frequency  non-selective  and  to  vary  slowly 
over  a  bit  period.  The  system  model  is  first  described,  followed  by  a  development  of  the 
error  probability  and  the  application  of  the  Gauss-Chebyshev  quadrature  rule. 

3.4.2  System  Description.  The  signal  set  for  this  NFSK  case  is 

s1(O  =  4,cos(ej1f  +  0)]  ^  ... 

1W  e  )  1  \[  0<t<Th.  (3.55) 

s2(t)  =  Accos(o)2t  +  0)\ 
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The  frequency  separation  between  oj\  and  a>2  is  assumed  to  be  large  enough  so  that  s\(i) 
and  52(0  occupy  non-overlapping  spectral  bands  [17],  Again,  it  is  assumed  that  the  two 
symbols  have  equal  a  priori  probability  and  the  AWGN  has  zero  mean  and  variance  NqH. 
The  Rayleigh  fading  is  again  assumed  to  be  frequency  non-selective  and  to  vary  slowly 
over  a  bit  interval.  The  receiver  used  for  this  model  is  shown  in  Figure  3.9. 


Figure  3.9.  NFSK  Receiver 


The  derivation  will  proceed  by  first  considering  the  case  where  there  is  no  Rayleigh 
fading  in  the  channel  and  then  proceed  to  account  for  the  fading.  The  two  possible 
received  signals  during  a  given  bit  interval,  as  defined  in  [17],  are 


*i(0 


— ^-cos  (a>xt  +  a)+ n(t) 

V  Tb 


0<t<Tb 


(3.56) 


J2E 

-j^cos(a>2t  +  a)+n(t) 


where  Eb  =  0.5 Ac2Tb,  a  is  the  random  phase  factor  that  is  uniformly  distributed  over 
(0,  2  tv),  and  n(t)  is  the  AWGN.  The  received  signal  set  is  then  resolved  into  the  space 
defined  by  the  orthogonal  basis  set 
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(3.57) 


With  the  basis  set  in  hand  and  assuming  that  si(t)  was  sent,  the  received  signal  is  defined 
by  the  data  vector  Z=  (X\,  X2,  Y\,  Y2)  where  the  components  are 

X\=  <jE^cosa  +  Nx„X2=Nx2 

I -  (y.JOj 

Y\  =  -yjEb  sin  a  +  Ny] ,  Y2=  Ny2 , 

where  Nx  1,2  and  Ny  \>2  are  uncorrelated  noise  components  with  zero  means  and  variances 
Ay 2.  Next,  the  PDFs  of  the  random  variables  R\  and  R2  defined  as 

R]  =VZ>2+7>2  (359) 


and 


K  =  V-rJ+n2  (3.60) 

are  sought.  Given  the  random  variable  a,  X\  and  Y\  are  Gaussian  random  variables  with 
means  of  (£’*)1/2cosa  and  -(Eb)m$\na,  respectively.  The  variance  of  both  X\  and  7i  are 
A/o/2.  Additionally,  X2  and  Y2  are  also  Gaussian  with  zero  means  and  variances  of  Ay 2. 
Because  X\  and  Y\  are  independent,  the  joint  PDF  of  X\  and  Fi  is  found  by  multiplying 
the  marginal  PDFs  of  X\  and  Y\  to  result  in 


1 


TT-ft,  -Jl^cosaJ  +(v,  +Je^ sin  at 
N0l 


(3.61) 


Similarly,  the  joint  PDF  of  X2  and  Y2,  which  does  not  depend  on  a,  is 


/V2(*2,72)  =  :^-exp 


(3.62) 
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The  PDFs  of  (3.61)  and  (3.62)  are  then  transformed  into  the  polar  coordinates  ri,2 
and  <j>\t2  defined  as 


x\,2  —  r\,2 1 


I— sin^, 
2  1 


1,2 


K  A 

y  i.2  ^y-^cos^ 


r,  2  >  0;  0  <  2  <  2;r 


(3.63) 


In  general,  a  variable  transformation  from  the  random  variables  U and  Finto  Z  =  g(U,V) 
and  W  =  h(U,V)  results  in  a  joint  PDF  in  terms  of  Z  and  W  defined  in  [15]  as 


Pzw 


(  \  PuAU l’Vl)  ,  PuAUn’Vn)  ,  , 

(z’w)=tt< — r+"‘T^7 — r+”'-  364 

\AuiM  V\un^A 


The  terms  un  and  v„  represent  the  real  roots  of  the  transformation  equations  and  J(u,v)  is 
the  Jacobian  of  the  transformation  defined  as 


j(w,v)  = 


In  this  transformation,  the  Jacobian  is 


dz 

dz 

du 

du 

du 

dv 

dz 

dw 

dw 

dw 

dv 

dv 

du 

dv 

dz 

dw 

(3.65) 


J(u,v)- 


—  sin  <j> 
2 


INn 


rcos</> 


hycostf 


N0  . 
-^-r  sin  <f>\ 


rNn 


(3.66) 


Applying  the  transformation  to  (3.61)  and  averaging  over  the  random  variable  a  results 


in 


(ri  >0\)  =  ~~  exP 


l(rd  +  2Eb) 

I0 

2  Eb  ) 
r„  — - 

1 

to 

0^ 

\"oJ 

r,  >  0,  0  <  (j>x  <  2n .  (3.67) 


Additionally,  the  joint  PDF  of  R2  and  <J>2  is 
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/V2(r2>^)  =  ”^exP 


r2>  0,  0  <  <j>2  <  In 


Finally,  integrating  (3.67)  and  (3.68)  over  <$1,2,  respectively,  results  in 
Pr,  fa )  =  r\  exP 


1 

f  2  2 Eb) 

\2Eb  1 

— 

h 

r,  J — - 

2 

l  ^0  J 

llNo) 

rx  >  0 


and 


PrM 2)  =  ^2  exp 


2 

2 


r2  >0. 


(3.68) 


(3.69) 


(3.70) 


From  (3.69)  and  (3.70),  R\  is  Rician  distributed  and  R2  is  Rayleigh  distributed.  A  similar 
derivation  in  the  case  where  S2(f)  was  sent  results  in  a  Rician  distribution  for  R2  and  a 
Rayleigh  distribution  for  R\.  Now,  with  the  unconditional  probabilities  of  Ri  and  R2  in 
hand,  the  error  probability  is  found. 


3.4.3  Probability  of  Error.  The  decision  about  which  signal  was  transmitted  is  based 
upon  the  random  variables  R\  and  R2.  Since  is  was  assumed  that  .s-i(/)  was  transmitted,  an 
error  occurs  if  R\  <  R2.  The  probability  that  R\  is  less  than  R2  is 


00  00 

p(r,  <R2\R])=  J pRi fa2 )dr2  =  jr2  exp 

Ri  Ri 

00 

=  |  exp(-  u)du  =  exp 


r  ..2  a 


v  2  j 


f  R^ 


\dr , 


V  2  J 


(3.71) 


where  a  variable  substitution  of  u  -  r2l2  was  made.  The  probability  in  (3.71)  is  then 
averaged  over  R\  to  obtain  a  probability  of  error  of 


pE = fp(p, 


_  e-E„/N 0 


jri  exp(-^2K 


u 


0 


v 


2  Eh 


N, 


I  drx . 


0  J 


(3.72) 
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Applying  the  integral 


'xe-ax2lAbx)dx  =  —ebllAa 
0V  ^  2  a 


from  [17]  to  (3.72)  obtains  the  closed-form  error  probability 


1 


(3.73) 


p  -  Pe-Eb/2N0 


(3.74) 


A  similar  derivation  in  the  case  where  ^(0  is  sent  results  in  the  same  error  probability  as 
in  (3.74);  therefore,  (3.74)  is  the  total  error  probability. 

The  next  step  in  the  derivation  is  to  consider  the  introduction  of  Rayleigh  fading 
to  the  channel.  The  inclusion  of  this  random  variable,  g,  to  the  model  results  in  a 
received  signal  of  the  form 


Zl,2  (0  ' 


1 


—  —  cos {yaX2t  +  a)+n(t)  0<t<Tb. 


(3.75) 


Since  g  is  a  Rayleigh  random  variable,  g1  will  have  a  chi-square  distribution  with  two- 
degrees  of  freedom.  Consequently,  the  signal-to-noise  ratio, 

Mg% 


Yb  =' 


Nn 


(3.76) 


will  also  have  a  chi-square  distribution  with  two-degrees  of  freedom  of  the  form 

p(yb)  =  -\-e  Yb 

Y\ 


(3.77) 


where  y’b  =  ybEig1)  is  again  the  average  signal-to-noise  ratio  as  in  the  BPSK  in  Rayleigh 
fading  case.  Finally,  averaging  over  yb  results  in  a  total  unconditional  probability  of  error 
of 
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oo  1  °° 

pb  =  \p(rb)PE\ydrb=—  jexP 


00 

f 

f  1 

1  Yl 

lexp 

-n 

— 1 - 

J 

0  ; 

r\)  j 

dYb 


1 


(3.78) 


2  Y\ 


1  1 

- 1 - — 

2  r 


2 +r\ 


b  J 


3.4.4  Application  of  the  Gauss-Chebyshev  Quadrature  Rule.  The  final  step  for  this 
case  is  to  again  find  the  two-sided  Laplace  transform  of  a  decision  metric.  First,  the 
transform  is  found  in  the  absence  of  Rayleigh  fading.  In  this  case,  the  decision  metric, 
assuming  that  R\  is  associated  with  the  signal  sent,  is 

Y  =  R2  -  R\ .  (3.79) 

As  shown  in  (3.79),  an  error  occurs  if  Y <  0,  which  corresponds  to  R2  >  R\2 .  Making  use 
of  the  PDFs  of  R\  and  Rj.  from  (3.69)  and  (3.70),  the  two-sided  Laplace  transform  is 

(5)  =  E^~sY  }=  E e ** } 

00  o°  i  /  \  r  i  ^3  80) 

=  J  f exp  [-  r,2  (0. 5  +  s)\0  (r,  ^ )it\  r2  exp[-r22(0.5-s)]s&-2 
0  Lo 

where  the  integration  could  be  performed  separately  because  Ri  and  R2  are  independent 
random  variables.  Applying  the  integral  of  (3.73)  to  the  inner  integral  of  (3.80)  results  in 


*,(?)=■ 


fi  1  11 

exp 

-n 

2(0.5 +s)j_ 

2(0.5 +  5) 

Next,  applying  the  integral  from  [18] 


j r2  exp[-  r22 (0.5  -  s)}/r2  .  (3.81) 


f  xme~axldx  = - , 

J  0/7(m+1)/2 


_r[(w+i)/2] 

2  a{ 


(3.82) 


where 


r(«+i)=»!, 


(3.83) 
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to  (3.81)  results  in 


<Ms)= 


fl  1  11 

exp 

-Yb 

^  2(0. 5  4" 

(3.84) 


4(0.5+5X0-5-^) 

Now  consider  the  case  where  the  Rayleigh  random  variable  g  again  affects  the  amplitude 
and  the  signal-to-noise  ratio.  Therefore,  the  result  in  (3.84)  is  conditioned  on  the  random 
SNR  Yb-  To  obtain  the  unconditional  Laplace  transform  requires  integrating  the  product 
of  the  conditional  transform  with  the  PDF  of  ybto  yield 


<D 


(s)  = 

—00  /  b  0 


eXP 


~Yb 


-  +  1- 


l  Y\  2(0.5  +  s)  JJ 


4(0.5  +  5X0.5  -5) 


dYb 


(3.85) 


Af„  (0.5  +  5X0.5  -  4-f  + 1  -  TZT"7~  \ 
{Yb  2(0.5 +  5)y 


The  error  probability,  using  (3.85)  and  the  Gauss-Chebyshev  quadrature  method, 
is  again  implemented  using  MATLAB  and  the  result  compared  to  those  obtained  from 
(3.78).  As  described  in  Chapter  II,  the  choice  for  the  constant  c  was  set  to  half  of  the 
smallest  positive  singularity.  Additionally,  the  error  probability  for  the  case  where  the 
channel  is  only  influenced  by  AWGN,  from  (3.74),  is  also  shown  for  reference.  The 
results  are  depicted  in  Figure  3.10.  As  shown  in  Figure  3.10,  only  ten  nodes  are  required 
to  obtain  a  maximum  percent  difference  of  less  than  one  percent  between  the  two 
methods.  The  fact  that  few  nodes  are  required  gives  comparable  computation  time 
between  the  two  methods.  The  code  is  shown  in  Section  4  of  Appendix  A. 
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Number  of  Nodes:  10;  Max.  Percent  Difference:  0.62% 


Figure  3.10.  Error  Probability  for  NFSK  in 
AWGN  and  Rayleigh  Fading 


CHAPTER  IV 


DISCUSSION  OF  THE  METHOD 

This  chapter  discusses  some  of  the  observations  made  when  using  the  Gauss- 
Chebyshev  quadrature  method  defined  in  Chapter  II  and  applied  in  Chapter  HI.  A  section 
discussing  some  general  observations  is  first  covered  followed  by  more  detail  on  how  to 
choose  the  constant  c. 

4.1  General  Observations 

The  utility  of  this  method  to  carry  out  error  probability  analysis  was  shown  in 
Chapter  III.  As  shown,  the  implementation  involves  defining  a  suitable  decision  metric 
and  then  finding  the  two-sided  Laplace  transform  of  the  PDF  associated  with  that  metric. 
The  method  used  to  obtain  this  transform  varied  slightly  between  the  different  models 
considered.  The  antipodal  baseband  signaling  case  involved  a  fairly  straightforward 
derivation  of  the  Laplace  transform.  Using  the  result  from  this  antipodal  case  in  the 
BPSK  in  AWGN  with  an  imperfect  carrier  reference  allowed  us  to  assume  a  Gaussian 
distribution  for  the  decision  metric.  Furthermore,  this  allowed  the  use  of  the  already 
known  form  for  the  Laplace  transform  of  a  Gaussian  PDF.  Then,  it  was  a  relatively  easy 
step  to  find  a  closed  form  Laplace  transform  by  integrating  the  product  of  the  conditional 
transform  with  the  phase  error  distribution.  Therefore,  it  is  seen  how  the  second  case 
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built  on  the  first.  A  similar  derivation  was  used  in  the  BPSK  in  AWGN  with  Rayleigh 
fading  case  but  a  closed  form  solution  was  not  obtained. 

Applying  the  technique  to  the  NFSK  case  where  the  channel  was  influenced  by 
AWGN  and  Rayleigh  fading  required  a  somewhat  different  approach.  A  closed  form 
solution  for  the  Laplace  transform  was  found  in  the  NFSK  case.  The  point  is  that  there  is 
no  single  way  of  easily  finding  the  Laplace  transform  of  the  PDF  associated  with  the 
decision  metric.  If  one  approach  toward  this  end  is  not  successful  in  obtaining  a  closed 
form  transform,  then  another  approach  may  be  in  order.  A  closed  form  solution  is 
obviously  desired  in  order  to  reduce  the  computational  complexity  when  implementing 
the  Gauss-Chebyshev  quadrature  numerical  technique.  Additionally,  the  utility  of  using 
results  from  one  case  to  solve  for  transform  of  another  case  shows  how  the  method  can 
save  time.  The  most  important  factor  in  obtaining  accurate  results  for  this  method 
revolves  around  finding  a  good  choice  for  the  constant  c,  which  is  now  discussed. 

4.2  Choosing  the  Constant  c 

The  most  important  factor  in  implementing  this  method  was  setting  the  constant  c. 
As  stated  in  the  method  development,  a  proper  choice  for  c  enables  accurate  results  to  be 
obtained  for  the  error  probability  while  using  the  fewest  nodes  possible.  This  obviously 
equates  to  requiring  less  computational  time.  To  recap,  it  was  stated  in  [2]  that  a  good 
choice  for  c  is  that  value  of  c  which  minimizes  <l>y(c).  This  corresponds  to  the  Chernoff 
bound.  Another  way  of  choosing  c,  as  stated  in  [3],  is  by  setting  it  to  half  the  value  of  the 
positive  singularity  of  OyfA)  which  is  closest  to  the  imaginary  axis.  In  most  cases,  setting 
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the  value  of  c  according  to  either  of  these  two  methods  will  result  in  accurate  results  with 
minimal  nodes. 

For  example,  consider  the  NFSK  in  AWGN  and  Rayleigh  fading  case.  The  two- 
sided  Laplace  transform  for  this  model  is 


from  (3.85),  repeated  here  for  illustration.  As  shown  in  (4.1),  the  positive  singularity  at 
s  =  0.5  is  that  which  is  positive  and  closest  to  the  imaginary  axis.  Therefore,  the  value  for 
c  in  this  case  was  set  to  0.25.  Indeed,  this  resulted  in  good  performance,  giving  a 
maximum  percent  difference  of  0.62%,  as  shown  in  Figure  3.10.  The  function  Oy(c)  for 
this  case  is  plotted  in  Figure  5.1  for  several  values  of  SNR.  This  value  of  c  worked  so 
well  because  at  c  =  0.25  the  function  <E>y(c)  is  well  behaved,  as  shown  in  Figure  4.1.  The 
integral  of  0}(.v)  is  thus  fairly  easy  to  approximate  using  the  Gauss-Chebyshev 
quadrature  rule  of  (2.38).  This  is  obvious  because  a  straight  line  is  much  easier  to 
approximate  than  say  a  quadratic.  For  comparison,  at  c  =  0.4,  <f>y(c)  is  less  well  behaved 
and  a  maximum  percentage  difference  of  6.7%  results  when  using  this  value  of  c  over  the 
same  SNR  range  as  when  c  =  0.25. 
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Figure  4.1.  <J> y{c)  for  NFSK  in  AWGN 
and  Rayleigh  Fading 


With  the  NFSK  example  in  hand,  it  can  be  stated  that  a  good  value  for  c  is  that 
value  for  which  <3>y(c)  is  well  behaved.  Essentially,  this  corresponds  to  where  the  slope 
of  3>y(c)  is  unchanging  or,  in  other  words,  where  the  second  derivative  of  Oy(c)  is  zero. 
This  hypothesis  is  tested  with  the  NFSK  case  because  it  is  evident  from  Figure  4.1  that 
c  =  0.25  may  not  be  a  good  choice  for  all  values  of  SNR.  Therefore,  a  technique  is 
implemented  to  obtain  values  of  c  where  the  second  derivative  of  ®y(c),  at  the  different 
SNR  values,  is  zero.  These  c  values  are  then  used  in  implementing  the  Gauss-Chebyshev 
quadrature  method.  There  was  no  appreciable  difference  between  results  obtained  when 
c  is  set  to  0.25  versus  c  chosen  as  described  here.  However,  when  the  c  optimization 
technique  is  applied  to  the  BPSK  in  AWGN  and  Rayleigh  fading  case,  performance  is 
improved  when  compared  with  other  c  values  tested. 
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In  the  end,  it  can  be  said  that  choosing  c  is  very  important  in  minimizing  the 
number  of  nodes  required  to  obtain  a  desired  degree  of  accuracy.  As  hypothesized,  a 
good  value  is  where  the  second  derivative  of  Oy(c)  is  closest  to  zero.  This  could  be 
potentially  found  in  closed  form  or  numerically.  However,  numerical  determination  will 
also  require  computing  time.  Therefore,  if  the  singularities  of  Oy(.s’)  are  easily 
identifiable,  then  choosing  c  to  be  half  the  value  of  the  real,  positive  singularity  closest  to 
the  imaginary  axis  seems  to  work  well.  If  a  singularity  does  not  exist  or  is  not  easily 
identifiable,  finding  the  second  derivative  of  Oy(c)  and  determining  where  it  is  zero 
results  in  a  good  choice  for  c.  The  next  best  technique  for  finding  good  c  values  is  by 
plotting  <1 >y(c)  and  graphically  attempting  to  determine  where  the  function  is  well 
behaved  to  choose  c.  This  worked  pretty  well  for  the  BPSK  case  in  AWGN  and  with  a 
phase  error  in  the  receiver.  However,  if  such  a  value  is  strongly  dependent  on  the  SNR,  a 
good  value  of  c  will  be  different  for  each  SNR  value  and  then  numerically  finding  good 
values  for  c  may  be  the  easiest  method.  However,  as  previously  stated,  this  will  add  to 
the  computational  time  required  to  implement  the  method. 


CHAPTER  V 


CONCLUSIONS 


5.1  Conclusions 

The  technique  described  in  [2]  to  analyze  error  probability  in  communication 
systems  using  the  two-sided  Laplace  transform  and  a  Gauss-Quadrature  rule  has  been 
fully  described.  Implementation  of  several  test  cases  has  verified  that  the  technique 
works  for  models  of  varying  complexity.  Methods  for  selecting  the  constant  c  have  been 
explained  and  tested.  These  methods  have  been  shown  to  provide  accurate  results  when 
compared  to  analytically  derived  error  probabilities  for  the  cases  considered. 
Additionally,  the  technique  was  implemented  assuming  equal  a  priori  probabilities  for 
each  binary  signal  set.  The  method  can  still  be  used  in  the  case  where  the  a  priori 
probabilities  are  not  equal.  Furthermore,  this  would  not  change  the  resulting  Laplace 
transforms. 

The  most  significant  utilization  of  this  technique  applies  to  models  where  closed 
form  expressions  for  the  error  probability  are  not  known.  Therefore,  in  many  cases, 
numerical  integration  is  required  to  obtain  the  error  probability.  As  shown  for  the  case  of 
BPSK  in  AWGN  with  phase  error  in  the  receiver,  the  error  probability  obtained  using  this 
method  resulted  in  significant  time  savings  over  numerically  integrating  the  traditional 
error  probability  expression  of  (3.30)  while  obtaining  similar  results.  Therefore,  it  is 
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suggested  that  when  a  closed  form  error  probability  can  not  be  found  analytically,  this 
method  may  allow  error  probabilities  to  be  computed  faster. 

5.2  Suggestions  for  Further  Research 

An  avenue  of  further  research  would  be  to  apply  this  technique  to  cases  where 
Oyfs)  is  not  analytically  known.  As  suggested  in  [2],  the  method  can  still  be  used  if  only 
a  few  numerical  values  of  Oy(s)  are  known.  Investigating  this  has  the  potential  to  allow 
analysis  of  error  performance  if  a  closed  form  expression  for  <Dy(s)  is  not  obtainable. 

Also,  further  application  of  this  technique  to  models  of  higher  complexity  shows 
good  potential.  For  example,  one  suggestion  is  to  apply  this  method  in  analyzing 
diversity  reception  of  DPSK  and  NFSK.  In  [22],  this  method  was  successfully  applied  to 
analysis  of  the  error  performance  of  wideband  direct-sequence  spread-spectrum  (DSSS). 

Investigating  the  possible  extension  of  this  technique  to  models  with  M- ary  signal 
sets  may  also  prove  to  be  useful.  Specifically,  using  this  method  to  analyze  QPSK  should 
be  a  relatively  easy  step.  The  effect  of  phase  error  at  the  receiver  could  be  considered  in 
a  similar  manner  as  that  in  the  BPSK  of  this  thesis.  Additionally,  one  could  consider 
possible  cross-talk  problems  that  arise  in  the  system.  The  method  could  be  further 
extended  toM-ary  PSK  and  other  M-ary  systems. 

Finally,  investigating  the  recursive  implementation  of  this  method  would  possibly 
reduce  computational  requirements.  Recursive  implementation  seems  possible  if  the 
nodes  are  doubled  in  each  successive  implementation.  This  is  because  the  abscissa 
locations  when  doubling  the  nodes  would  be  maintained.  Additionally,  implementing  the 
method  when  there  is  no  known  error  probability  could  be  conducted.  In  such  a  situation, 
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the  method  could  be  recursively  implemented  until  the  error  probability  change  between 
implementations  falls  below  some  set  threshold.  Another  way  to  ensure  the  method  is 
working  would  be  to  compare  the  results  to  those  obtained  through  simulation. 
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Appendix  A 

MATLAB  PROGRAMS 

This  Appendix  shows  the  code  used  to  generate  the  various  error  probabilities  for 
the  different  models  explored  in  this  paper.  Commands  used  to  plot  the  results  are  not 
displayed.  Note  that,  with  the  exception  of  the  antipodal  case,  the  programs  are 
implemented  with  loops.  Further  computational  time  would  be  saved  if  these  loops  were 
implemented  in  matrix  form  like  for  the  antipodal  case,  but  that  was  not  the  focus  of  my 
efforts. 


A.l  Antipodal  Baseband  Signaling 


The  MATLAB  program  used  to  find  the  error  probability  in  the  antipodal 
baseband  signaling  model  is  considered  here.  The  program  is  given  below: 

clear; 

A=2;  T_b=.5;  %Amp.  of  Signal;  Bit  Time 

n=30;  m_Y=A*T_b;  %#  of  Nodes  to  be  Used ;  m_Y\+A  Sent 
k=l:n/2;  %Index  for  tau__k 

tau_k=tan ( (2*k-l) *pi/n) ; 

SNR__dB=l :  15 ; 

SNR=10. A (SNR_dB/10) ; 

var_Y=AA2*T_bA2 . / (2*SNR) ;  %Variance  of  Y 
c=(AA2*T_bA2) ./ (2*var_Y) ;  %Setting  for  c;  equal  to  SNR 
PHI=exp ( ( ( .5*var_Y.*c. A2)  1  * ( 1+j *tau_k)  .A2)- 
m_Y*c ' * ( 1+j  *tau_k) ) ; 
x=ones ( 1, length ( SNR) ) ; 

%Creates  a  Matrix  of  tau__k  (s)  to  avoid  using  a  loop  to 
per  from  the  summation  for  several  SNR  values 
tau_k= ( tau_k 1 *x)  *; 

PE_GQ= ( 1/n) *sum ( ( real ( PHI ) +tau_k . *imag ( PHI ) ) 1 ) ; 

%P_E  using  Gauss-Quad .  Method 
PE ( SNR_dB) = . 5*erf c ( sqrt { SNR) ) ; 

%Known  Closed  Form  Error  Probability 
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A.2  BPSK  in  AWGN  with  Tikinov  Phase  Error 


The  program  for  finding  the  error  probabilities  for  the  BPSK  case  is  shown  here. 
The  following  part  of  the  program  shows  general  variables  and  constants  used  in  both 

methods  used: 

clear ; 

A=2 ;  Tb=.5;  %Amp .  of  Signal;  Bit  time 

n=100;  m_Y=A*Tb ;  %#  of  Nodes;  m__Y\+A  Sent 
k=l:n/2;  %Index  for  tau_k 

tau_k=tan ( (2*k-l) *pi/n) ; 
track=-10 ; 

%The  amount  of  difference  desired  between  tracking  loop 
SNR  (in  dB) and  AWGN  SNR .  (Positive  values  mean  tracking 
loop  SNR  is  greater  than  the  AWGN  SNR) 

SNRdB=l:15;  SNR=1 0 . A ( SNRdB . / 1 0 ) ; 

var_Y= (AA2*TbA2) . / (2*SNR) ;  %Variance  of  Y  given  phi 
rho=10 . A ( (SNRdB+track) . / 10 ) ; 

The  next  portion  of  the  program  shows  numerical  integration  using  MATLAB’s  quad8 
command. 

tic 

for  i=l : length (SNR) ; 

PE_Q8 (i) =quad8 ( 'petikinov* , -pi, pi, [ ] , [ ] , rho (i) , SNR (i) ) ; 
end; 

timel=toc; 

The  following  shows  the  function  definition  of  petikinov  called  above  in  the  quads 
command: 

function  y=petikinov (phi , rho, SNR) 

y= ( 1 . / ( 4*pi*besseli (0, rho ) ) ) . *exp ( rho . *cos (phi ) ) . * 
erf c (sqrt (SNR) ,*cos (phi) ) ; 

Finally,  the  Gauss-Quadrature  method  is  implemented  as  follows: 

tic 

for  i=l : length ( SNR) ; 
c=5 ; 

ww=c+j  *c. *tau_k; 

PHI= (besseli (0, (-A*Tb*ww+rho (i) ) ) . / 

(besseli  (0,  rho  (i )  )  )  )  .  *exp  (  .  5*var__Y  (i  )  *ww.  A2  )  ; 

PE_GQ ( i ) = ( 1/n) *sum ( real ( PHI ) +tau_k . *imag ( PHI ) ) ; 
end; 
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time2=toc; 

SAVE=round ( 100- ( time2/timel )  *100)  ; 

%Calculation  of  time  savings 
PE= . 5 . *erf  c ( sqrt ( SNR) ) ; 

%Prob.  of  Error  with  phase  error  =  0 


A.3  BPSK  in  AWGN  with  Rayleigh  Fading 


Again,  the  first  portion  of  the  program  involves  defining  constants  and  variables 

clear;  Ac=l;  Tb=l;  n=50;  k=l:n/2;  tau_k=tan ( (2*k-l) *pi/n) ; 

SNRdB=l : 15;  SNR=10 . A ( SNRdB/10 ) ;  N0=AcA2*Tb . / (2*SNR) ; 
c=linspace (.01,3,20) ; 

The  next  portion  of  the  program  involved  finding  a  good  value  for  the  constant  c  to 
obtain  accurate  results  with  the  fewest  possible  nodes.  This  involves  finding  that  portion 
of  the  function  Oy(c)  that  behaves  nicely.  This  enables  the  Gauss-Quadrature  method  to 
approximate  the  function  more  accurately  with  a  set  amount  of  nodes  or  allows  you  to 
decrease  the  number  of  nodes  used  to  any  desired  degree  of  accuracy  and  reduce 
computation  time.  The  code  is: 


tic; 

phic= [ ] ; 

for  jj=l : length (SNRdB) ; 
for  ii=l : length (c) ; 


s=c (ii) ; 
xx=sA2/2 ; 

yy=quad8 ( * fade’ , 0, 1, [ ] , [ ] , xx) ; 
uu=l-exp ( s A2/4 ) * ( s*sqrt (pi/4) - (sA2/2) *yy ) ; 
phic (length ( c) * ( j  j -1 ) +ii ) =exp ( . 5*1/ ( SNR ( j  j )  ) *Tb*s . A2 ) 
end; 
end; 

phi c=re shape (phic, length ( c) , length ( SNRdB) ) ; 

[ YY, opt] =min (abs (dif f (phic, 2 ) ) ) ; 

c_opt=c (opt) ; 

timel=toc; 


*uu; 


The  function  call  to  fade  in  the  quad8  command  calls  the  function  that  defines  the 
argument  of  the  numerical  integration  required  in  this  model  (reference  (3.54)).  The  code 


in  this  file  is: 


function  y=f ade ( z , w) 
y—exp (- . 5*w*z . A2 ) ; 
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Finally,  the  Gauss-Quadrature  method  is  applied  as  follows: 


tic; 

for  ii=l : length (NO ) ; 

^s=c+j  *  c*tau_k; 

s— c_opt  (ii)  H- j  * c _ opt  (ii)  *tau__k; 

xx=s . A2/2 ; 

for  j  j=l : length (xx) ; 

yy ( j  j ) =quad8 ( 1  fade ' ,  0 , 1 , [ ] , [ ] , xx ( j  j ) ) ; 
end; 

uu=s*sqrt (pi/4) - ( (s . A2) /2) . *yy; 

PHI=exp ( . 5*N0 (ii) *Tb*s . A2)  . *  (1-exp (s . A2/4)  . *uu) ; 

PE_GQ (ii ) = ( 1/n) *sum ( real ( PHI ) +tau_k . *imag ( PHI ) ) ; 
clear  s; 
end; 

time2=toc; 

PE= . 5* erf c (sqrt (SNR) ) ; 

PE_T=. 5* ( 1-sqrt (SNR. / ( 1+SNR) ) ) ;  ^Theoretical,  closed  form  error 
probability 

DIFFERENCE= ( 1/ 100 ) * round (max ( abs ( PE_GQ- 
PE_T) ./ ( .5* (PE_GQ+PE_T) ) *100) *100) ; 

%Calculates  the  percent  difference 


A.4  NFSK  in  AWGN  with  Rayleigh  Fading 


Finally,  the  NFSK  program  is  presented.  The  first  portion  of  the  program  again  deals 
with  declaring  the  various  constants  and  variables  used  in  the  program. 

clear; 

n=10;  t;#  of  Nodes 

k-l:n/2;  %  Index  for  tau__k 

tau_k=tan ( (2*k-l) *pi/n) ;  ^Location  of  abscissas 
SNRdB=l : 20 ; 

SNR=10. A (SNRdB/10)  ; 

Next,  the  theoretical  error  probability  in  the  case  of  an  AWGN  channel  (PE)  and  the 
close  form  solution  for  the  AWGN  channel  with  Rayleigh  fading  (PE  FADEl)  are 
computed. 


PE=. 5*exp (- . 5*SNR) ; 

PE  FADE1=. 5* ( 1 . / ( 1+ . 5*SNR) ) ; 
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The  Gauss-Quadrature  method  is  now  applied  with  the  constant  c  set  to  half  of  the  value 
of  the  smallest  positive  singularity. 


for  i=l : length (SNR) ; 

c=.25;  %Sets  the  optimum  Value  of  c 

ww=c+j*c. *tau  k; 

A= ( 1/ ( SNR ( i ) )T-(1./ (2*(.5+ww) ) )+l; 

B=4*SNR(i) * ( .5-ww) .* ( . 5+ww) ; 

PHI=1 . / (A. *B) ; 

PE_GQ (i) = ( 1/n) * sum (real ( PHI ) +tau_k . *imag ( PHI ) ) ; 
end; 

ERROR= 1 0  0  *  ab  s ( ( P E_FADE 1 - P E_GQ ) )  . / ( . 5* ( PE_FADE1+PE_GQ) ) ; 
MAX_ERROR=round(100*max (ERROR) ) /100; 

MAX  ERROR=num2str (MAX  ERROR) ; 
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Next,  consider  the  second  integral  of  (B.3).  This  could  easily  be  represented  in  the  form 
of  a  complementary  error  function  but,  since  the  limit  is  complex,  MATLAB  could  not 
implement  it.  Therefore,  the  second  integral  is  put  in  the  form 


Applying  the  integral 

from  [18]  to  the  first  integral  in  (B.5)  and  making  a  change  of  variables  of  x  =  (a/b)z  to 
the  second  integral  results  in 


Finally,  the  original  integral  is  expressed  in  the  form 


(B.8) 


